Spin- wave analysis of the transverse-field Ising model on the checkerboard lattice 



(N 

o 

(N 



i 

-a 

o 
o 



> 

O 

(N 



X 



Louis-Paul Henry, 1 Peter C. W. Holdsworth, 1 Frederic Mila, 2 and Tommaso Roscilde 1 

' Laboratoire de Physique, CNRS UMR 5672, Ecole Normale Superieure de Lyon, 
Universite de Lyon, 46 Allee d'ltalie, Lyon, F-69364, France 
2 Institut de Theorie des Phenomenes Physique, Ecole Polytechnique Federate de Lausanne (EPFLj, CH-1015 Lausanne, Switzerland 

(Dated: February 8, 2012) 

The ground state properties of the S = 1/2 transverse-field Ising model on the checkerboard lattice are 
studied using linear spin wave theory. We consider the general case of different couplings between nearest 
neighbors (Ji) and next-to-nearest neighbors (J2). In zero field the system displays a large degeneracy of the 
ground state, which is exponential in the system size (for Ji = J2) or in the system's linear dimensions (for 
Ji > Ji). Quantum fluctuations induced by a transverse field are found to be unable to lift this degeneracy 
in favor of a classically ordered state at the harmonic level. This remarkable fact suggests that a quantum- 
disordered ground state can be instead promoted when non-linear fluctuations are accounted for, in agreement 
with existing results for the isotropic case Ji = J2. Moreover spin-wave theory shows sizable regions of 
instability which are further candidates for quantum-disordered behavior. 

PACS numbers: 75.30.Ds, 75.30.Kz, 75.50.Ee, 75.10.Jm 



I. INTRODUCTION 

Frustrated quantum magnets represent one of the richest 
playgrounds to investigate quantum collective phenomena 1 . 
Indeed known models of frustrated quantum magnets might 
admit quantum ground states without any classical counter- 
part, such as valence bond crystals, or resonating-valence- 
bond spin liquids 2 . Most of the investigations have focused 
on Heisenberg antiferromagnets, due to their relevance to real 
compounds, although the lack of well controlled analytical 
or numerical approaches for the bulk properties of these sys- 
tems still leaves the question open on the true nature of their 
ground state 3 . On the other hand, quantum-disordered ground 
states have been shown to emerge in models with strongly 
anisotropic interactions (in spin space and/or in real space) 
which appear to provide the first controlled realizations of 
quantum spin liquids with a topological natural—. Some of 
these systems are related to quantum dimer models^, which 
have provided the first known examples of quantum spin- 
liquid ground states^. 

In this class of anisotropic systems a special role is played 
by frustrated quantum Ising models, namely Ising systems en- 
riched with quantum fluctuations, coming either from an ex- 
change coupling in the transverse spin components, or from 
a transverse magnetic fiel d 11 ' 12 . In the absence of quantum 
fluctuations, Ising models on frustrated lattices have generally 
a classical spin liquid nature, namely they exhibit an expo- 
nential degeneracy of the ground state - as it is the case for 
the Ising antiferromagnet on the triangular lattice, kagome lat- 
tice, checkerboard lattice, pyrochlore lattice, etc., with how- 
ever correlations that can be long-ranged, algebraic or short- 
ranged. The effect of quantum fluctuations is in general that of 
lifting the large degeneracy of ground states, leading either to 
the emergence of an ordered ground state (as e.g. in the case 
of the triangular lattice-^) or to a ground state with novel spin- 
liquid properties (as it appears to be the case for the pyrochlore 
lattice^). In this study we focus on the Ising antiferromagnet 
on the checkerboard lattice, which represents a fundamental 
model of frustrated magnetism in two dimensions. Indeed, 



in the case of spatially isotropic interactions the ground state 
properties of the system are equivalent to those of the unbi- 
ased six-vertex model, also known as square icei^. The effect 
of quantum fluctuations on such a system has been the sub- 
ject of several recent investigation a 1 J i 14 ~ .1£, focusing particu- 
larly on the limit of weak quantum fluctuations, treated within 
degenerate perturbation theory. 

In the present study we adopt a different strategy, which 
allows us to treat arbitrarily strong quantum fluctuations in a 
generalized version of square ice. We investigate the 5=1/2 
Ising model on the checkerboard lattice with different cou- 
plings along the coordinate axes ( J{) and along the diagonals 
(J2), as shown in Fig.Q] The anisotropy in the couplings al- 
lows to introduces a bias in the vertex weights of the corre- 
sponding vertex model, reducing the degeneracy to exponen- 
tial in the linear dimensions of the system ( J 2 > J\), or even 
to a finite value (for J\ > J 2). The application of a transverse 
field allows therefore to investigate the effect of quantum fluc- 
tuations on a classical manifold of states with variable degen- 
eracy. Quantum fluctuations are investigated via linear spin- 
wave theory, accounting for the harmonic fluctuations around 
the classical ground state. Albeit limited to the harmonic ap- 
proximation, such an approach allows for investigation of ar- 
bitrarily strong fields, and its breakdown signals the candidate 
regions in the phase diagram where novel quantum disordered 
behavior can be expected. On the side J2 > J\ the infinite 
degeneracy of the ground state would a priori make the spin- 
wave analysis impossible, given the exceedingly large number 
of possible classical reference states. In fact we demonstrate 
that the spin-wave spectrum does not depend on the partic- 
ular classical ground state chosen as a reference. This im- 
plies that the spectrum of linear excitations is defined unam- 
biguously in the J2 > J\ region - this situation persists also 
in the isotropic case of square ice, J\ = J2, for which the 
lowest branch of the excitation spectrum is a flat band. But 
this implies as well that the classical degeneracy remains un- 
altered in presence of harmonic quantum fluctuations, which 
means that only non-linear quantum fluctuations can lift the 
degeneracy, a situation previously encountered in a number of 
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FIG. 1: Checkerboard lattice. Full lines represent couplings with 
strength Ji , dashed lines represent couplings with strength J2 . 



frustrated magnets^—. Moreover, harmonic quantum fluc- 
tuations triggered by a sufficiently strong field are found to 
reverse the classical hierarchy of ordered states close to the 
isotropic (Ji = J 2) limit, suggesting that anharmonic fluc- 
tuations might destabilize the classical order completely. Fi- 
nally the classical order parameter is found to be completely 
washed out by quantum fluctuations for strong fields close to 
the classical polarization transition, and for couplings close to 
the isotropic case. We can therefore conclude that the trans- 
verse field Ising model (TFIM) on the checkerboard lattice can 
harbor quantum disordered phases for strong frustration and 
fields, and that non-linear quantum fluctuations are expected 
to play a major role in the case of extensive degeneracy of the 
ground state. 

The paper is structured as follows: Sec. [TT] describes the 
model, its behavior in zero applied field, and the classical be- 
havior in a transverse field; Sec.[lII]reviews known results on 
quantum square ice from degenerate perturbation theory and 
numerics; Sec. |IV] describes linear spin-wave theory as ap- 
plied to the various regimes of the checkerboard lattice TFIM; 
Sec.[V]discusses the phase diagram emerging from spin-wave 
theory; and finally conclusions are drawn in Sec. I VII along 
with a discussion about physical realizations. The general 
framework of spin-wave theory for the TFIM on arbitrary lat- 
tices is presented in Appendix |Al while the spin-wave ob- 
servables for the checkerboard lattice are discussed in Ap- 
pendix [B] 



II. CLASSICAL BEHAVIOR OF THE TRAS VERSE-FIELD 
ISING MODEL ON THE CHECKERBOARD LATTICE 



A. Model Hamiltonian and ground-state properties in zero 
field 



The Hamiltonian of the TFIM on the checkerboard lattice 
reads 

u = ,h ]T s? s| + ,h J2 sfs* - r ]T sf . (i) 



where Si are quantum spins of length S, satisfying \Si\ 2 = 
S(S + 1) and [Sf , Sf ] = ie af)l S]. The first sum in Eq. Q 
runs over the nearest-neighbor bonds of a square lattice, while 
the second sum runs over the next-to-nearest-neighbor (diag- 
onal) bonds on a staggered array of plaquettes (see Fig. Q]). 
We consider here frustrated antiferromagnetic couplings J\, 
J2 > 0. T is a transverse magnetic field, introducing quan- 
tum fluctuations in the system. As we will see in the fol- 
lowing section, in zero field the ground-state properties of the 
above model are equivalent to those of an m-vertex model 
with m = 2, 4 or 6 depending on the Hamiltonian parame- 
ters. Motivated by this equivalence, in the following we will 
indicate as vertices (denoted by M) the squares with additional 
diagonal J2 couplings, and as plaquettes (denoted by □) the 
squares without diagonal couplings. 

When r = we can easily rewrite the Hamiltonian in the 
following form 



TL — J2 Ziicc + {Ji — J2) Si<Sj 



if J 1 > J2, and 



% — J\ h[ CC + (J2 — Ji) SfS| 



(2) 



(3) 



if J 2 > J 1, where we have introduced the square-ice Hamil- 
tonian 



^-icc ^ ^ 



(4) 



With this choice the ground state is identified in two steps: 1) 
firstly, one needs to impose on each spin the constraint that 
Sf = ±S, and on each vertex the zero-magnetization con- 
straint, = X^ieKi &i = ; these two constraints mini- 
mize the first term on the right hand side of both Eq. (O and 
(01. The Mb = constraint corresponds to the so-called ice 
rule for square ice, and therefore we will hereafter denote the 
states which satisfy it (and which satisfy Sf = ±S) as ice- 
rule states; 2) secondly, one needs to choose, among the ice- 
rule states, those corresponding to a minimum of the second 
term in the r.h.s. of Eqs. (f2| and (0); this term is always anti- 
ferromagnetic by construction, so that it can be minimized by 
zero-magnetization vertices. 

In the special case J\ = J2 we recover % = J\(2)hi ce , 
namely the square ice model, whose ground states are only 
constrained by the ice rule (along with the condition Sf = 
±S). The ground state properties of this system are equiva- 
lent to the 6-vertex model, displaying an exponential degener- 
acy of the ground stated. This degeneracy is (partially) lifted 
when Ji 7^ J2 . 

If Ji > J2 we need to select among all zero-magnetization 
vertices with satisfied antiferromagnetic J\ bonds. This fur- 
ther constraint imposes that each vertex must take the antifer- 
romagnetic (Neel) configuration in Fig. 0a) or its spin-flipped 
version (2-vertex model), and it reduces the possible ground 
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FIG. 2: Crossed plaquettes that obey the "ice rule", (a) On the 
Neel plaquettes, the nearest neighbor links are satisfied, (b) On the 
"collinear" plaquettes, the next nearest neighbor links are satisfied. 



states to the two Neel-ordered states. The ground-state energy 
per spin takes the value i?Neei = — (2Ji — J2)S 2 . 

If on the contrary J 2 > Ji, the zero-magnetization vertices 
must satisfy the J 2 bonds, and therefore they take one of the 
4 collinear configurations in Fig. |2|a), reducing the ground 
state of the system to that of a 4-vertex model. The collinear 
ground states (namely containing only collinear vertices) are 
massively degenerate, because flipping a linear chain of spins 
along the J 2 diagonals does not alter any of the constraints 
which the ground state must satisfy. Therefore the ground- 
state degeneracy is equal to 2 Nd , where Nd is the number of 
J 2 diagonals. The ground-state energy per spin takes the value 

-Ecoll = —J2S . 

In summary the degeneracy of the ice-rule states is partially 
lifted when J% 7^ J 2 , and ice-rule states are organized into a 
band with energy width (per spin) given by 

w = |-En&i — E co u\ = 2|Ji — J2\S 2 . (5) 

The above discussion is valid for an arbitrary value of the 
spin length S, and in particular also in the classical limit S — > 
00, in which one can introduce continuous spins Si = Si/S 
with unit length \Si\ 2 = 1. On the other hand, we will see 
that the application of a transverse field introduces significant 
differences for the ground-state properties depending on the 
spin length S. 



B. Excited states in zero field 

The nature of the excited states is dependent upon the spin 
length. In the following we will concentrate on the S = 1/2 
case, which will be the main focus of this paper. In this limit, 
when Ji > J 2 the lowest energy excitations correspond either 
to: 

1) a single spin flip, costing an energy 

A = (2 - n) J x + (1 - vi) J 2 (6) 

where v p is the number of frustrated J p links (p = 1,2) con- 
nected to each site. For Neel states, V\ = and v 2 = 2, while 
for collinear states v\=1 and v 2 = 0. Notice that for generic 



ice-rule states one has v\ + v 2 = 2. 

2) a joint flip of all the spins on a plaquette (plaquette flip). 
This operation has the lowest energy when applied to aflip- 
pable plaquette, with the property that its flip connects the 
initial ice-rule state to another ice-rule state. All the neigh- 
boring vertices of a flippable plaquette share with the pla- 
quette a bond with zero magnetization, so that the plaque- 
tte flip will not alter the vertex magnetization. This imposes 
that the flippable plaquette has a local Neel configuration. In 
Neel states all plaquettes are flippable, with an energy cost of 
Api a q = 4(Ji — J2). In the collinear states, only a portion 
of the plaquettes are flippable (at most one half as in the state 
depicted in Fig. 2]i, with an energy cost of Api aq = 4(J 2 — Ji). 

Comparing the energy cost of a plaquette flip with that of 
a spin flip, we find that plaquette flips are the lowest-energy 
excitations in the parameter range 2/3 < J 2 j J\ < 4/3, and 
outside of this range spin flips are instead the excitations with 
the lowest energy. 



C. S — > 00 limit in a transverse magnetic field 

In this section we discuss the effect of a transverse field 
on classical continuous spins (S — > 00). The correspond- 
ing ground state configurations will serve as a template for 
the spin-wave analysis in the quantum case. The transverse 
field introduces a canting of the spins along the x axis by 
an angle #. This angle increases with the intensity of the 
field up to a critical value T c at which $ = ir/2, corre- 
sponding to the polarization of the spins along x. For both 
the collinear and the Neel ground states, each spin sees the 
same local field (in modulus) created by the neighboring spins 
(h ioc = ±(AJ 1 -2J 2 )Sz for Neel states, andh ioc = ±2J 2 Sz 
for collinear states). Therefore the application of the external 
field will create the same canting angle i? on each spin. Upon 
canting, the spin configuration becomes Sf = S sin ■& and 
S* = iScostf. 




d 



FIG. 3: Classical rotation angle induced by the transverse field. 

The classical energy per spin admits the compact expres- 
sion, valid for both the Neel and collinear states: 

e d = [{vi -2),h + {v 2 -l)J 2 ]S 2 cos 2 tf-rSsintf 
= [u 2 {J 2 - Ji)- J 2 ]S 2 cos 2 d-TSsind 

(7) 

Minimizing the energy per spin with respect to we find 

si ntf = min(— T — -^l) (8) 

\2S [J 2 - v 2 ( J 2 - Ji)\ J 

If T > T c = 25 [J 2 — v 2 (J 2 — Jx)], the system becomes 
completely polarized in the transverse direction. For T > T c 
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FIG. 4: Reference classical ground states. The circles indicate flip- 
pable plaquettes. 



the classical ground-state energy per spin takes the value 

e cl = -S 2 [J 2 - u 2 ( J 2 - Ji)] (1 + sin 2 0) . (9) 

The resulting classical phase diagram for S — > oo is shown in 
Fig.H 
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FIG. 5: Classical phase diagram of the transverse-field Ising model 
on the checkerboard lattice. 



tive constant): 
UcS = ~ A3 51 

flippablc □ 
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5fc 5; 



(10) 

where (ijkl) are the four sites on a plaquette (in clockwise or- 
der), and the sum runs on flippable plaquettes. Various claims 
exist in the literature 1 14 that this Hamiltonian will lift the de- 
generacy of ice-rule states in favor of the Neel state, based 
on the fact that the Neel state has the largest number of flip- 
pable plaquettes. In fact, this Hamiltonian is completely off- 
diagonal for the ice-rule states, so that it cannot favor a spe- 
cific ice-rule state, but only a resonant superposition thereof. 
In particular a flippable plaquette, once flipped, turns its four 
neighboring vertices from Neel configurations to collinear 
configurations or viceversa, namely the effective Hamiltonian 
will resonantly connect local Neel configurations with local 
collinear ones. At the same time, two corner-sharing flippable 
plaquettes cannot be both flipped without leaving the ice-rule 
manifold, which suggests that Neel-collinear resonances will 
be localized to single plaquettes. A rough approximation to 
the ground state of Eq. (IToT > is therefore a state in which 
a checkerboard subset of non-corner sharing plaquettes res- 
onate between two flippable states, giving rise to a resonating- 
plaquette solid (RPS) 

l*o> w I*rps> = I] (I UjtkU) + 1 4.it*4*ti» /V5 
□ 

(ii) 

where the primed product runs on a sublattice of plaquettes. 
This state breaks the two-fold symmetry between the two pla- 
quette sublattices. While direct numerical investigations of 
the square ice model in a transverse field are not known to 
us, there exists in the literature a series of numerical studies 
of square ice with different quantum perturbations, which all 
map perturbatively onto the effective Hamiltonian of Eq. ( fTOb - 
A direct simulation of Eq. (fT0T > is reported in Ref. |T5| while 
Refs. [lrjl Fnl and[l8l focus on hardcore bosons on the checker- 
board lattice at half filling, with strong nearest-neighborrepul- 
sion, and weak hopping (with either positive or negative sign) 
between nearest neighbors and next-to-nearest neighbors. The 
latter model is equivalent to square ice with weak ferromag- 
netic/antiferromagnetic exchange terms for the x and y spin 
components. The common result of all these studies is that 
the ground state for weak quantum perturbations has indeed 
long-range RPS order, and no magnetic order. Therefore we 
expect spin-wave theory to break down or become inconclu- 
sive in this limit - which is indeed one of the main results of 
the following analysis. 



III. SQUARE ICE IN A TRANSVERSE FIELD : RESULTS 
FROM PERTURBATION THEORY 



IV. LINEAR SPIN- WAVE THEORY 



For square ice (Ji = J 2 ), in the case T <C Ji, Ji the trans- 
verse field can be treated via degenerate perturbation theory 
on the manifold of ice-rule states. Considering all terms up 
to fourth order—, one can easily find the following effective 
Hamiltonian for the subspace of ice-rule states (up to an addi- 



In the following we describe a treatment of quantum fluctu- 
ations introduced by the transverse field in the 5=1/2 case 
based on a linear spin wave expansio n 25,27 . We will then treat 
separately the spectrum of excitations above the various clas- 
sical reference states of the system: Neel, collinear, and fully 
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polarized. 

A. Spin-boson transformation 

We begin by considering a generic classical ground state 
with long-range magnetic order, and with a magnetic unit cell 
containing n spins. We denote Si ; „ the p-th spin (p = 1 . . . n) 
of the Z-th cell. As seen in Section HTCl in the classical limit an 
applied transverse field rotates the spins around the y-axis by 
an angle We introduce a local rotation of the spin configu- 
ration, Si.p = <T p TZy((Tp , d)Si t p, where a p = 1( — 1) if the spin 
in zero field has positive (negative) projection along the z axis, 
and lZ y (±i?) is the rotation matrix of an angle ±$ around the 
y axis. In the classical limit 5 — >• oo the ground state is a sim- 
ple ferromagnetic state for the Si tP spins, namely Sf p = S 
everywhere. 

We then consider small quantum fluctuations around this 
classical reference state, by transforming the quantum spins to 
bosons via a linearised Holstein-Primakoff transformation 28 
valid in the limit of a small number of bosons n^ p -C 25: 

Sf, P = s ~ a \. P a hP S?, P ~ \j\ (4, P + a i,p) ■ ( 12 ) 



Here ai lP and aj are bosonic operators, satisfying 

k P .4p] = land [ a S> a S] = a 

B. Harmonic Hamiltonian for ordered ice-rule states 

The Hamiltonian is then expanded up to quadratic order in 
the bosonic operators (the linear terms vanish by construc- 
tion). In the following we will specialize the discussion to 
reference classical states which in zero field are ice-rule states 
with long-range magnetic order, namely the states which min- 
imize the energy in the classical limit 5 — > oo, and whose 
ordered structure allows to build a spin-wave theory. These 
states have the property that the number of frustrated bonds of 
type 1 and 2, v\ and v-x, is the same for every site. Under these 
generic assumptions the quadratic bosonic Hamiltonian reads 

HlSW = A^cl + JlUnn + J-lUnnn + TH r (13) 

with 



H nn = 25 v 2 cos 2 1? ^ a l, P a i,p + 2 sin2 ^ X] \ a l,p a i',p' + a \,p a i',p' + hx -) 
i,p (ip,l'p') 

Hnnn = 25(1 - I> 2 ) COS 2 § ^ a l,p a l,P + ^ ^ ^ \ a l,p a l',p' + a l, P a l'-P' + hx -) (14) 

l,p i(lp-l'p')) 

H r = -Tsm§y^ j a] v ai iP 



This is a remarkable result, in that the spin-wave Hamilto- 
nian depends only on the frustration parameters v\ and v%, 
while it is completely independent of the geometry of the 
unit cell. The frustration parameters v\ , V2 distinguish among 
Neel states and collinear states, but they are not able to distin- 
guish among different collinear states. At the square-ice point 
Ji = J2 the dependence on the frustration parameter drops. 
Therefore, as we will discuss further, quantum corrections at 
the harmonic level are not able to lift the degeneracy between 
ordered ice-rule states, regardless of the size of their magnetic 
unit cell. This result can be extended even to disordered ice- 
rule states, which can be regarded as ordered ones with an 
infinite unit cell. 

To diagonalize the spin-wave Hamiltonian, Eq. (TT4T >. we 
first introduce the Fourier transform of the bosonic operators 
and then perform a ?i-modes Bogoliubov transformation as 
described in Appendix lAl The Hamiltonian then becomes 

U = N (e cl u: Kp U p b k , p + (15) 

k,p ^ ' 



where e = 25 [J2 — ^2 {J2 — Ji)} ■ 



C. Neel state 

Let us first consider the Neel state, defined for T < 
25(2Ji — J2). Its unit cell contains n = 2 spins (see fig- 
ure HJ. The unit cells form a rotated square lattice with vec- 
tors x' = (x — y) and y' — (x + y). The diagonalization 
of Mk shows that the spectrum of the magnon excitations is 
gapped whenever the classical Neel state is defined - i.e. if 
r < 25(2 Ji — J2). Its lower band is plotted in figure|6] It has 
minima at (0, 0) and at the four corners of the Brillouin zone. 
This corresponds to the structure of the classical Neel state. 



D. Collinear states 

As already mentioned in Sec. IIV Al all collinear states ad- 
mit the same frustration parameter U2, and hence the same 
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FIG. 6: Lowest band of the magnon spectrum for T — Ji/2 and 
various values of J2/J1, around the Neel state. For the purpose of 
readability, the bands have been offset by the energy of their lower 
edge. The lower band edge (corresponding to the minimum excita- 
tion gap) is plotted in the lower panel. 
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FIG. 7: Lowest band of the magnon spectrum for T = Ji/2 and var- 
ious values of J2/J1, around any collinear state. For the purpose of 
readability, the bands have been offset by the energy of their lower 
edge. The lower band edge (corresponding to the minimum excita- 
tion gap) is plotted in the lower panel. 



spin-wave Hamiltonian. This means that they possess the 
same spectrum of harmonic spin-wave excitations (but folded 
into a smaller Brillouin zone, the larger the unit cell), and 
that zero-point quantum fluctuations cannot lift the degener- 
acy among them. 

We will then specify the discussion to the particular 
collinear state represented in figure H] Its unit cell contains 
n = 8 spins; the unit cells form a square lattice with vec- 
tors x' = 2 (x — y) and y' = 2 (x + y). While not being the 
simplest of all collinear states, this state is relevant because 
it can be energetically stabilized against other collinear states 
by e.g. dipolar interactions, which are relevant for realistic ice 
models^. The magnon dispersion relation, obtained by diago- 
nalizing Mjj, is shown in Fig. [7] It shows a finite gap, and two 
lines of minimum-energy degenerate modes along the axes of 
the first Brillouin zone of the magnetic lattice (1/8 of the Bril- 
louin zone of the geometric lattice, shown in Fig. |7). These 
degenerate modes traveling with momentum (k x ,±k x ) for all 
k x values can be associated with deconfined monopole-like 30 
pairs, obtained by flipping a finite string of spins along a J2- 
diagonal of the checkerboard lattice. Given the degeneracy 
of all collinear states, not perturbed by quantum fluctuations, 
these pairs are deconfined along the J 2 -diagonals, and their 
energy is independent of momentum as long as it satisfies the 
constraint of diagonal motion. 



E. Quantum square ice 

At the square-ice point J\ = J2 = J, spin-wave theory 
built around any ice-rule state with canting spins produces 
the same excitation spectrum and zero-point quantum fluc- 
tuations. Therefore harmonic quantum fluctuations are not 
able to lift the degeneracy of the classical ice-rule manifold, 
and the elementary excitations are identical to the classical 
case, namely deconfined monopole pairs moving with arbi- 
trary momentum. This is reflected in the spin-wave disper- 
sion, showing a perfectly flat band for the lowest energy ex- 
citations, and with a gap equal to the classical value, namely 
A cl = 2e = 2JS 2 . 



F. Polarized states 

For large T, the classical reference state is totally polar- 
ized along the field. Nonetheless quantum fluctuations intro- 
duce deviations from the polarized state, given that the field 
term in the Hamiltonian does not commute with the Ising cou- 
plings. We build a spin-wave Hamiltonian around the polar- 
ized state analogously to what has been done for the ice-rule 
states. Even though in the classical spin configuration all the 
spins have the same orientation, the magnetic unit cell con- 
tains two sites (exactly as in the case of the Neel state), due to 
the fact that the checkerboard lattice is not a Bravais lattice. 
The bosonic excitations correspond to deviations of the spin 
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J 2 /J 1= 0.7 -r= 1.43 Jg/J =0.7-T= 1.3 




FIG. 8: Upper panels: lowest band of the magnon spectrum around 
the polarized state for T > T c (left column) and for r = T c (right 
column). Lower panels: Minimum excitation gap: (left) as a function 
of Ja/Ji for T/Ji = 1.3 (green line) and for T/Ji = 1.43 (blue 
line); (right) as a function of T / J\ for J2/J1 = 0.7 and 1.3 (the 
curves for both cases coincide). 



from full polarization along the x axis. 

The magnon spectrum is shown in Fig. [8] It displays softer 
modes at the four corners of the square-lattice Brillouin zone 
(for Ji > J 2 ), and along the edges of the checkerboard-lattice 
Brillouin zone (for J\ < J 2 ). These modes become gapless 
when approaching the critical field T c , signaling the instability 
of the fully polarized state to a Neel state (for J\ > J 2 ) and to 
degenerate collinear states (for J\ < J 2 ). 



V. RESULTS OF THE SPIN WAVE ANALYSIS 

The main observables from linear spin-wave theory are rep- 
resented by the internal energy E = (Hlsw) and the order 
parameter m = (S z )/S, whose expressions are given in Ap- 
pendix|B] These two quantities allow us to extract the quantum 
phase diagram of our system in the harmonic approximation. 

For each value of J 2 and I\ we can a priori choose be- 
tween three families of reference states as candidate ground 
state: the Neel states, the fully polarized state and the degen- 
erate collinear states. Notice that classically, the Neel and the 
collinear states differ in the value of the correlations between 
next-to-nearest neighbors C' 2 ' = (5fS , |)^ ij -^ = (yi — 1). 

For the Neel state, = S 2 , whereas for all collinear states 

^collinear = — ■ ^ natural definition of the ground state (or 
ground-state manifold) identified by spin wave theory is the 



state (or the family of states) which has the lowest energy, 
whose order parameter is finite and whose next-to-nearest- 
neighbor correlations have the correct sign; moreover the sta- 
bility of the state requires also that the spin-wave frequencies 
be real numbers. The satisfaction of these four conditions al- 
lows us to identify the phases indicated in Figs. [9] and [10] 
which correspond to the classical phases for the same param- 
eter ranges. 

A clear region of instability of spin-wave theory is found 
close to the classical transition line between the fully polar- 
ized phase and the Neel and collinear phases, as indicated by 
the black region in Fig. [10] In this parameter range the or- 
der parameters for all reference states are found to vanish, as 
shown in Fig. [12] clearly signaling the onset of a quantum- 
disordered phase. 




0.5 



°0 05 1 1.5 2 

FIG. 9: Phase diagram of the 5=1/2 transverse-field Ising model 
on a checkerboard lattice from spin-wave theory. The hatched, and 
dense-hatched regions correspond to the next-to-nearest neighbor 
correlation inversion region and to the energy-hierarchy inversion re- 
gion, respectively (see text). Below the dotted lines one can find - 
in the classical limit - two states which are (local) energy minima 
and which have order in the z spin components: one state with Neel 
order, and the other with collinear order. Above this line one of the 
two states becomes the polarized state. Therefore collinear and Neel 
ordered states can only be compared energetically below the dotted 
line. 



However, a number of other anomalies revealed by the 
linear-spin wave theory close to the transition to the param- 
agnetic phase point to a significantly larger region where the 
classical behaviour is probably destroyed by quantum fluctu- 
ations. 

First of all, a qualitative deviation from the classical behav- 
ior is observed when approaching the square ice limit J% = J 2 
in a strong field 0.7 < T < J-y. In this range (indicated by the 
dense-hatched region in Fig. [9]) we observe an inversion of 
the energy hierarchy between Neel and collinear states with 
respect to the classical case. This occurs despite the fact that 
quantum fluctuations are stronger for the energetically favored 
phase, as shown by the order parameters of the two phases 
obeying an opposite hierarchy (namely m^ e i > 77i co m near 
when Ersieei > -^collinear, and viceversa, see Fig. [TTJ. The 
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inversion in the energy hierarchy is due to strong quantum 
corrections to the classical energy, which change qualitatively 
the dependence of both Neel and collinear energies on J2 / J\ . 

This strong quantum effect of energy hierarchy inversion 
suggests that classical order might be unstable around the hi- 
erarchy inversion region in Fig. [9] when considering quantum 
fluctuations beyond linear spin-wave theory. The real ground 
state of the system may then be an intermediate phase which 
cannot be described within the linear spin wave approxima- 
tion. 

Another strong quantum effect is revealed close to the clas- 
sical phase boundaries. While classically E'n&i and -Ecoiiinear 
are monotonic functions of J2/J1, they become non- 
monotonic around the above mentioned field range. In par- 
ticular i^Neei grows with increasing J2 / Ji until it reaches a 
maximum, beyond which it starts to decrease; from a classi- 
cal point of view this is quite surprising. According to the 
Hellmann-Feynman theorem, the next nearest neighbor cor- 
relations C^ 2 ' are given by the derivative of the energy with 
respect to J2, namely C' 2 -* = d(H)/dJ2- Consequently, a 
change of sign in the derivative of E corresponds to a change 
of sign in which means that the harmonic ground state 
is dramatically different from the reference state. The locus 
of the maxima in the energy as a function of J2 / Ji represents 
the lower bound of the hatched region in Fig. [9] and Fig. [10] 
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FIG. 10: Zoom of Fig.[9]around the classical tri-critical point. In the 
black domain, the order parameters of all considered phases vanish. 
Other marked zones are as in Fig. [9] 



A further element of inconsistency of spin-wave theory is 
offered by the apparent violation of the Hellman-Feynman 
theorem 



d(H) 
dT 



(16) 



The transverse magnetization, defined as (m x ) = (S x ) is 
proportional to m in the harmonic approximation : (S x ) = 
S sin m. The deviation from Eq. (TToT l comes along with a 
strong non-monotonic behavior, as shown in Fig. [13] This be- 
havior is unphysical, implying a negative susceptibility. This 
signals again that an ordered reference state does not lead to 




1 




0.8 - 
0.6 - 
0.4 - 
0.2 - 

S U - 3 J 2 /J, 

FIG. 1 1 : Average energy (upper panel) and order parameter (lower 

panel) associated with the Neel (red solid line) and collinear (blue 
solid line) reference states. The dashed lines correspond to the clas- 
sical energies. Both panels are for V / Ji — 0.9. 



consistent results. We identify therefore an additional region 
of strong anharmonic fluctuations with the magnetization dip 
in the m x vs. T curve. This dip lies inside the region [Ti, T2] 
din 

where (T±) = andm x (T2) = m x (Ti). The latter con- 
dition marks the upper bound of the hatched region in Figs. [9] 
andUO] 

Indeed Eq. ([ToT i ceases to be strictly valid in the linear spin- 
wave approximation because in general the linear spin-wave 
Hamiltonian does not have the form H = Hq — TS X . From 
inspection of Eqs. dT3l l and ( fT4b it is immediate to see that each 
term of the Hamiltonian depends on V through 6(F), except 
for the polarized case in which 9 = tt/2 independently of the 
value of r. In fact, to recover the identity of Eq. ( [ToT i at order 
1/S, one has to include corrections to the harmonic ground 
stated. Even though Eq. ( [Tol l ceases to be an identity within 
linear spin-wave theory, the precision with which Eq. (IT6b is 
approximately verified can be used as a further criterion for 
the validity of the linear spin-wave approximation. In Fig [13] 
(S x ) is compared to d(H) /dT: upon increasing the field, the 
quantities show a significant deviation from each other when 
approaching the region already identified before as showing 
significant inconsistencies of spin-wave theory. 

A strong violation of the Hellman-Feynman theorem 
Eq. ( fT6l is observed also in the case J2 = 0. In this limit 
the system reduces to the square-lattice Ising model in a trans- 
verse field, which features a well-known quantum phase tran- 
sition (for a field T c / Ji w 1.5- see Ref.l32l) between Neel or- 
der and a quantum paramagnetic state. Interestingly the devi- 
ation between m x and d(H) / dT builds up when approaching 
the quantum-critical field T c , showing that spin-wave theory 
is able to signal the quantum phase transition via the break- 
down of the consistency of its results with known theorems. 
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FIG. 12: Order parameter vs field for different values of J2. It van- 
ishes at the critical value of the field and a small gap opens for J2 
close to Ji where none of the three states has a finite order parame- 
ter. The insets are closer views around the critical value of the field. 
The hatched and dense-hatched regions correspond those identified 
in the phase diagram (Fig.[9j(. 




FIG. 13: Field-induced magnetization as a function of the field for 
different values of J2/J1. The highlighted regions are as in Fig.l 121 
The dotted lines correspond to the opposite of the derivatives of the 
average energies with respect to the magnetic field T. These curves 
deviate from the average induced magnetization, signaling a viola- 
tion of the Hellmann-Feynman theorem (see text). 



the next-to-nearest-neighbor coupling J2 exceeds the nearest- 
neighbor one J 1, and which is exponential in the system size 
for J\ = J2, corresponding to the square-ice limit. This im- 
plies that spin wave theory is inconclusive regarding the ques- 
tion of which classical ground state is selected by quantum ef- 
fects, and that non-linear quantum fluctuations play a central 
role in lifting the degeneracy in the exact ground state. This 
result is consistent with existing studies of the square ice in a 
weak transverse field, for which degenerate perturbation the- 
ory and numerics suggest a quantum-disordered ground state 
with the structure of a resonating plaquette solid. 

Our results suggest that a quantum-disordered ground state 
persists beyond this limit. Spin-wave theory indicates an 
anomalous inversion in the classical hierarchy between refer- 
ence states for strong fields and for J\ ~ J2, suggesting that 
quantum fluctuations beyond the harmonic approximation can 
destroy classical order in the system in that parameter range. 
Moreover spin-wave theory breaks down completely close to 
the classical transition line to full polarization, where linear 
quantum fluctuations are able to suppress all spin components. 
This identifies an interesting candidate region for quantum- 
disordered ground states. 

Recent experimental developments have led to controlled 
realizations of frustrated Ising models, both in the classical 
limit and in the presence of tunable transverse fields. In par- 
ticular the classical square-ice model, enriched with long- 
range dipolar interactions, is realized in recent experiments on 
square-lattice arrays of nano-patterned magnetic domains 33 . 
While quantum tunneling of the magnetization in these sys- 
tems is not realistic, given the mesoscopic size of the magnetic 
moments, one can envisage scaling down the components 
of these artificial ice systems to single-molecule magnets 34 , 
which can be arranged into regular arrays on a surface 3 -^. 
Transverse-field Ising models are currently realized by arrays 
of trapped ions^ 3 ^, where the internal states of the ions can 
be coupled with Ising Hamiltonians via virtual phonons, and 
where transverse fields are created by Raman laser schemes. 
The ions can be individually trapped by micro-traps, which 
can in turn be arranged into arbitrary planar arrays 3 ^, encom- 
passing the checkerboard geometry explored here as well as 
other frustrated structures. Therefore the theoretical investi- 
gation of transverse-field Ising models on frustrated lattices is 
very compelling, as it promises to lead to the realization of 
novel quantum states in controlled artificial spin systems in 
the near future. 
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In this study we have applied a linear spin-wave theory 
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is exponential in the linear dimension of the system when 
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Appendix A: Spin- wave theory for general Ising Hamiltonians 
in a transverse field 

We present here several general formulas to study the 
quadratic quantum fluctuation in a generic transverse field 
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Ising system. We consider a generic classical ground state 
with long-range magnetic order, and with a magnetic unit cell 
containing n spins. We denote Si iP the p-th spin (p = 1 . . ,ri) 
of the /-th cell. The most general Hamiltonian supporting such 
a ground state has the form 



H 



TFI 



\pp' ^l,p^V ,p' 



ip,i f p' 



(Al) 

Here r\ is the position of a reference site in the I— th unit cell, 
and J(Ar) is a n x n matrix containing the couplings between 
spins in unit cells at a distance Ar. 

In the classical limit an applied transverse field rotates the 
p-th spins around the y-axis by an angle t) p . As already dis- 
cussed in Section IIV Al we introduce a local rotation of the 
spin configuration, Si. p = a p TZ y (a p 'dp)Si t p, where a v = ±1 
is the orientation of the spin in zero field. The rotation has the 
effect of reducing the S = oo ground state to a perfectly fer- 
romagnetic one. The classical energy of the p-th spin of each 
cell has the expression 

S 2 

£ci, P = — o-pcos^p ^2 [J(Ar)] pp , oy cos $ p > -ST smtf p , 

Ar,p' 

so that the total classical energy can be written as 

Ed = — ^ £cz.p = — Tre ci (A2) 



where we have introduced the matrix [e c i]p,p' = s c i,pS PtP '. 

We then consider small quantum fluctuations around this 
classical reference state, by transforming the quantum spins to 
bosons via a linearised Holstein-Primakoff transformation 28 
valid in the limit of a small number of bosons m tP <C 2S 1 : 



(A3) 



Here ai iP and a\ p are bosonic operators, satisfying 

[aj lP , a\ ] — 1 and [aj^, a^ p ] = 0. The angles i? p are chosen 
so that the classical reference state is stable. Thus the linear 
terms in the bosonic operators vanish. 

The quadratic Hamiltonian then reads 



H 2 = E d +^2 hp a i P a ip 
i, P 



(A4) 



\ Yl J( ri ' - r 'W ( a l + ai p) ( a «V + a i'p' 



where 



ip,i ; p f 



hp = 2e c i iP /S + r sinz?p 

J(Ar) P pi = J(Ar) P pi sin -dp sin $ p / 



(A5) 



ators and of the interaction 



a k.j 



J(k) = ^c- ikA "J(Ar) 



(A6) 



The quadratic Hamitonian can then be written in the compact 
form 

H 2 = ^ £ - t) + \ £ (AT) 



where 



Zip = — e c i tP + T sin i?p 

A k = ( fl kl.---i a kn' a -k.lr--,1-M) 

A k A k \ _ fe cl n \ 
A k A ky l ^0„ e ch 

-- ifj(k) + J(k)t 



M k 
A 



This Hamiltonian can be diagonalized by a n- 
mode Bogolyubov transformation. This consists 

in finding the transformation yl k = T k B k , with 
B k = (6 kl ,...,fe kn ,fo_ k a,...,6_ kj „) T , such that 

A{M k A k = ^2tul p) bl p b Kp and [& k , P ,& k , P ] = 1 and 



We introduce the matrix S, given by 

y _ fin On 
\0n -It. 

the matrix Z k of the right eigenvectors of XM k , and 
the unitary matrix U k such that U^Z^lZ k U k = 

diag(£ k , . . . , l k ) = L k . The transformation matrix 
T k is then obtained as^i^ 2 - 



T k = Z k U k \L k \-^ 



(A9) 



(v) 

In particular, the eigenmodes w k are the eigenvalues of 
SM k . 



If the matrices A k and e c i commute (which is the case for 
the Neel and collinear states of the checkerboard Ising model 

(v) 

studied here, having e c i = £oI n ), the eigenmodes Lu k can be 
expressed in terms of the eigenvalues A k of A k in the form 



,0) 



We then introduce the Fourier transform of the bosonic oper- 




(A10) 
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Appendix B: Energy and magnetization in the case of Neel and 
collinear states 

As mentioned above, in the particular case of the Neel and 
collinear states of the checkerboard lattice Ising model studied 
in this work, the classical energies e c i. p are all equal to Eq. The 
mean energy of the system then reads 



(E) = N(£ d -£ /2) + -^2uj ktP 



(Bl) 



As Ak is proportional to T 2 , so are its eigenvalues. We can 
then expand the mean energy per spin in powers of T 2 . We 
will introduce rescaled eigenvalues Ak. p defined as Ak. p = 
£oT 2 Ak,p. 



2/// 



k rn—1 



-^a m Tr ( — 



k.rn 



£0 



(B2) 



where the a m are defined by y/1 + Ax — 1 = 4 a m x m . 

ra—l 

In the case of a Neel or a collinear case, the trace of Ak av- 
erages to zero in the Brillouin zone. Thus the first non-zero 



correction to the classical energy in Eq. ( IB2I ) is of fourth order 
inT. 

In both cases, if we introduce the ratio J = J2/J1, we get 
(e)=e ct +Y,c p (J)( -, Jl ^ Yt 2 p (B3) 



p>i 



S{V2 ~ Jf 



where the c p coefficients only depend on J, and not on the 
considered state. It is then obvious that the expansion be- 
comes independent of the classical state if J = 1. Thus the 
classical degeneracy of the ice model is not lifted by harmonic 
fluctuations. 

Similarly, if all the Wk, P are real (which is the case whenever 
spin-wave theory holds), we have 



m = 1 — 



^-E(4 (P) 4 P) 



NS 



k.p 



( 2Wk !f : 



4NS e Q 



— y 



2uj 



k.p 



(B4) 



1 x — > x — > -~ — ' m n 



m— 1 k,i 



where the /? m are defined by \/l + 4x 

00 
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